Photoproduction of Charmonia and 
Total Charmonium-Proton Cross Sections 



J. Hiifner a ' 6 , Yu.P. Ivanov a ' 6 ' c , B.Z. Kopeliovich 6 ' c and A.V.Tarasov a ' 6,c 



a 



Institut fur Theoretische Physik der Universitdt, Philosophenweg 19, 69120 Heidelberg, Germany 
k Max-Planck Institut fur Kernphysik, Postfach 103980, 69029 Heidelberg, Germany 
Joint Institute for Nuclear Research, Dubna, 14-1980 Moscow Region, Russia 

February 1, 2008 

Abstract 

Elastic virtual photoproduction cross sections 7*p — > J/tp(ip')p and total charmonium- 
nucleon cross sections for J/tp, tp 1 and x states are calculated in a parameter free way 
with the light-cone dipole formalism and the same input: factorization in impact pa- 
rameters, light-cone wave functions for the 7* and the charmonia, and the universal 
phenomenological dipole cross section which is fitted to other data. The charmonium 
wave functions are calculated with four known realistic potentials, and two models for 
the dipole cross section are tested. Very good agreement with data for the cross section 
of charmonium photoproduction is found in a wide range of s and Q 2 . The inclusion 
of the Melosh spin rotation increases the ip' photoproduction rate by a factor 2 — 3 and 
removes previously observed discrepancies in the ip' to J ftp ratio in photoproduction. 
We also calculate the charmonium-proton cross sections whose absolute values and 
energy dependences are found to correlate strongly with the sizes of the states. 
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1 Introduction 



The dynamics of production and interaction of charmonia has drawn attention since their 
discovery back in 1973. As these heavy mesons have a small size it has been expected 
that hadronic cross sections may be calculated relying on perturbative QCD. The study of 
charmonium production became even more intense after charmonium suppression had been 
suggested as a probe for the creation and interaction of quark-gluon plasma in relativistic 
heavy ion collisions [|TJ. 

Since we will never have direct experimental information on charmonium-nucleon total 
cross sections one has to extract it from other data for example from elastic photoproduction 
of charmonia 7p — > J/ip(ip') p. The widespread believe that one can rely on the vector 
dominance model (VDM) is based on previous experience the with photoproduction of p 
mesons. However, even a dispersion approach shows that this is quite a risky way, because 
the Jfij) pole in the complex Q 2 plane is nearly 20 times farther away from the physical region 
than the p pole. The multichannel analysis performed in || demonstrates that the corrections 
are huge, cx^f p turns out to be more that three times larger than the VDM prediction. 
Unfortunately, more exact predictions of the multichannel approach, especially for if)', need 
knowledge of many diagonal and off-diagonal amplitudes which are easily summed only if 
one uses the oversimplified oscillator wave functions and a gg-proton cross section of the 
form (J q q(rT) oc r^, where tt is the transverse qq separation. 

Instead, one may switch to the quark basis, which should be equivalent to the hadronic 
basis because of completeness. In this representation the procedure of extracting cr^f p from 
photoproduction data cannot be realized directly, but has to be replaced by a different strat- 
egy. Namely, as soon as one has expressions for the wave functions of charmonia and the 
universal dipole cross section a q q(rT,s), one can predict both, the experimentally known 
charmonium photoproduction cross sections and the unknown cx^f^ p . If the photopro- 
duction data are well described one may have some confidence in the predictions for the 
Gtlot^ ' >P - Of course this procedure will be model dependent, but we believe that this is the 
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best use of photo-production data one can presently make. This program was performed for 
the first time in || . The aim of this paper is not to propose a conceptually new scheme, but 
to calculate within a given approach as accurately as possible and without any free param- 
eters. Wherever there is room for arbitrariness, like forms for the color dipole cross section 
and those for for charmonium wave functions, we use and compare other author's proposals, 
which have been tested on data different from those used here. 

In the light-cone dipole approach the two processes, photoproduction and charmonium- 
nucleon elastic scattering look as shown in Fig. [I] || . The corresponding expressions for the 




Figure 1: Schematic representation of the amplitudes for the reactions 7*p — > ipp 
(left) and ipp elastic scattering (right) in the rest frame of the proton. The cc 
fluctuation of the photon and the ip with transverse separation tt and cm. energy 
-y/i interact with the target proton via the cross section a(rx, s) and produce a 
J/ip or ip' . 



forward amplitudes read 

1 

M rp (s,Q 2 ) = £ (da (d 2 f T ^\a,r T )a qq (r T ,s)^\a,f T ,Q 2 ) , (1) 

1 

M+ P {s) = £ fdafd 2 f T ^\a,r T )a q ,(r T ,s)^\a,r T ). (2) 
w o 

Here the summation runs over spin indexes fj,, p, of the c and c quarks, Q 2 is the photon 
virtuality, $ 7 *(a, ry, Q 2 ) is the light-cone distribution function of the photon for a cc fluc- 
tuation of separation tt and relative fraction a of the photon light-cone momentum carried 
by c or c. Correspondingly, <3>,/,(a, rr) is the light-cone wave function of J/ip, ip' and x ( om y 
in Eq. 0). The dipole cross section (J qg -(r T , s) mediates the transition (c/ Fig. p. 
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In Section |2] we review the status of the factorized light-cone approach to photoproduction 
of heavy quarkonia. Besides the well known distribution function of quarks in the photon, it 
needs knowledge of the universal flavor independent dipole cross section which depends on 
the transverse qq separation and energy. In Section [2.1| we introduce two parameterizations 
available in the literature. 

Making use of the nonrelativistic approximation for heavy quarkonia in Section 172 we 
solve the Schrodinger equation with four types of relativistic potentials available in the 
literature. The next most difficult step is a Lorentz boost to the infinite momentum frame 



discussed in Section p73| . Although this procedure is ill defined and no unambiguous recipe is 
known, we apply the standard and widely used one. We put a special emphasis on importance 
of the Melosh spin transformation, which turns out to be very important. 



The final expression for the photoproduction cross sections is presented in Section 3.1 



and results are compared with available data for J/ip production in Section \5.2[ Although 
the calculations are parameter free they demonstrate a very good agreement with data. 
The ratio of ip' to J / ip photoproduction yields has drawn attention recently since previous 



calculations grossly underestimate experimental values. It is demonstrated in Section 3.3 



that the Melosh spin transformation which has been overlooked previously, and accompanies 
the Lorentz boost may be the reason. It has a dramatic impact on the ip' photoproduction 
increasing its yield by a factor 2 — 3 in a good agreement with data. 

After we will have demonstrated that the approach under discussion quantitatively ex- 
plains the photoproduction data, we calculate in Section |] the total charmonium-nucleon 
cross sections for J/i[),i/;' and x~ s - We predict quite a steep energy dependence for these 
cross sections slightly varying for different charmonia. Although the cross sections correlate 
with the mean charmonium size, this dependence is slower than cx (r|>), and this fact finds 
a simple explanation. In Section ^ we compare our estimates for charmonium-nucleon cross 
sections with the effective absorption cross section of charmonium which can be extracted 
from data on nuclear attenuation of J/ip and ip' . Agreement is rather good. 
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Our results are summarized in Section |6] where we also discuss the physics of energy 
dependence of the cross sections and the status of our approach. Special attention is given 
to nuclear attenuation of charmonia which is affected by formation and coherence time 
phenomena in an important way. 

2 Light-cone dipole formalism for virtual photopro- 
duction of charmonia off nucleons 

The light cone variable describing longitudinal motion which is invariant to Lorentz boosts 
is the fraction a = p^ /pi* of the photon light-cone momentum pi = E^* + p 7 * carried by 
the quark or antiquark. In the nonrelativistic approximation (assuming no relative motion 
of c and c) a = 1/2 (e.g. ||), otherwise one should integrate over a (see Eq. ([]])). For 
transversely (T) and longitudinally (L) polarized photons the perturbative photon-quark 
distribution function in Eq. ([I]) reads (|, ||] , 

^f(a,r T ,Q 2 ) = ^f^Z cX fd T , L ^K {er T ) , (3) 

I 7T 

where 

Xc = iv y X*c ; ( 4 ) 
X and x are the spinors of the c-quark and antiquark respectively; Z c = 2/3. Kq^tt) is the 
modified Bessel function with 

e 2 = a(l - a)Q 2 + m 2 c . (5) 
The operators Ot,l have the form: 

T = m c a ■ e 7 + i(l — 2a) (a ■ n) (e 7 • V rT ) + (n x e 7 ) ■ V rT , (6) 

6 L = 2Qa(l-a)a-n , (7) 

where n = p/p is a unit vector parallel to the photon momentum and e is the polarization 
vector of the photon. Effects of the non-perturbative interaction within the qq fluctuation 
are negligible for the heavy charmed quarks. 
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The color dipole cross section o qq - (tt, s) is poorly known from first principles. It is 
expected to vanish oc r\ at small rx —* due to color screening and to level off at 
large separations due to a finite range of gluon propagation. We employ phenomenological 
approaches described in section p7T| . 

The charmonium wave function is well defined in the rest frame where one can rely on 
the Schrodinger equation. We present solutions for four potentials proposed in the literature 
(section |2.2|) . As soon as the rest frame wave function is known, one may be tempted to 
apply the Lorentz transformation to the cc pair as it would be a classical system and boost it 
to the infinite momentum frame. However, quantum effects are important and in the infinite 
momentum frame a series of different Fock states emerges from the Lorentz boost. (Compare 
with a Lorentz boost of a positronium: Weizsaker- Williams photons appear.) Therefore the 
lowest |cc) component in the infinite momentum frame does not represent the |cc) in the 
rest frame. We rely on the widely used procedure for the generation of the light-cone wave 



functions of charmonia and describe it in section 1273. 



2.1 Phenomenological dipole cross section 

The dipole formalism for hadronic interactions introduced in expands the hadronic cross 
section over the eigen states of the interaction which in QCD are the dipoles with a definite 
transverse separation (see (0)). Correspondingly, the values of the dipole cross section a qq -{r T ) 
for different are the eigenvalues of the elastic amplitude operator. This cross section is 
flavor invariant, due to universality of the QCD coupling, and vanishes like (J qq (rT) oc r\ for 
tt—i-O. The latter property is sometimes referred to as color transparency. 

The total cross sections for all hadrons and (virtual) photons are known to rise with en- 
ergy. Apparently, the energy dependence cannot originate from the hadronic wave functions 
in Eqs. ([I], 0), but only from the dipole cross section. In the approximation of two-gluon 
exchange used in || the dipole cross section is constant, the energy dependence originates 
from higher order corrections related to gluon radiation. On the other way, one can stay 
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with two-gluon exchange, but involve higher Fock states which contain gluons in addition 
to the qq. Both approaches correspond to the same set of Feynman graphs. We prefer to 
introduce energy dependence into <J q q(rT, s) and not include higher Fock states into the wave 
functions. 

For small size dipoles essential for DIS one may apply perturbative QCD and the energy 
dependence comes as an effect of of gluon radiation treated in the leading-log(l/x) approx- 
imation |7|, §]. In the opposite limit of large separations typical for light hadrons one can 
also calculate the effects of gluon bremsstrahlung making use of smallness of the quark-gluon 
correlation radius 0. 

However, the intermediate case we are interested in, is the most complicated one as 
usual. No reliable way to sum up higher order corrections is known so far. Therefore we use 
a phenomenological form which interpolates between the two limiting cases of small and large 
separations. Few parameterizations are available in the literature, we choose two of them 
which are simple, but quite successful in describing data and denote them by the initials of 



the authors as "GBW" JTCj and "KST" JTT 
We have 

"GBW": a g g(r T ,x) = 23.03 



1 _ e -r%/rl{x) 



mb , 



x \ - 144 
r (x) = 0.4 ( — ] fm 



where Xq = 3.04 ■ 10~ 4 . The proton structure function calculated with this parameterization 
fits very well all available data at small x and for a wide range of Q 2 [pTD| . However, it 
obviously fails describing the hadronic total cross sections, since it never exceeds the value 
23.03 mb. The rr-dependence guarantees Bjorken scaling for DIS at high Q 2 , however, Bjorken 
x is not a well defined quantity in the soft limit. Instead we use the prescription of ||12| , 
x = (Ml + Q 2 )/s, where is the charmonium mass. 



This problem as well as the difficulty with the definition of x have been fixed in [|1T 
The dipole cross section is treated as a function of the cm. energy a/s, rather than x, since 
■s/s is more appropriate for hadronic processes. A similarly simple form for the dipole cross 
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section is used 



'KST": 



(9) 



The values and energy dependence of hadronic cross sections is guaranteed by the choice of 



cr [s) 
ro(s) 



23.6 ( — 

.So 



0.08 



1 + 



3r2(s 



<ri); mb ' 



0.88 



so 



-0.14 



fm . 



(10) 
(11) 



The energy dependent radius fo(s) is fitted to data for the proton structure function F$(x, Q 2 ), 
So = 1000 GeV 2 and the mean square of the pion charge radius (r 2 h ) = 0.44 fm 2 . The im- 
provement at large separations leads to a somewhat worse description of the proton struc- 
ture function at large Q 2 . Apparently, the cross section dependent on energy, rather than 
x, cannot provide Bjorken scaling. Indeed, parameterization (|9]) is successful only up to 
Q 2 « 10 GeV 2 . 

In fact, the cases we are interested in, charmonium production and interaction, are just 
in between the regions where either of these parameterization is successful. Therefore, we 
suppose that the difference between predictions using Eq. @ and @ is a measure of the 
theoretical uncertainty which fortunately turns out to be rather small. 

We demontrate in Fig. || a few examples of ^-dependence of the dipole cross section 
at different energies for both parameterization. The KST cross section reveals a nontrivial 
behavior, it rises with energy at tt < 3 fm 2 , but decreases at larger separations. This is 
however a temporary effect, <7o(s) reaches minimum at y's « 77 GeV and then slowly rises 
at higher energies. Such a peculiar behavior is a consequence of our original intention to 
reproduce the energy dependence of the hadronic cross sections a t g t oc s 08 keeping the form 
(|) of the cross section. Of course data are insensitive to the cross section at such large 
separations. 

Both GBW and KST cross section vanish oc r\ at small ry, however considerably deviate 
from this simple behavior at large separations. Quite often, the simplest parameterization 
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Figure 2: The dipole cross section as function of r\ at energies y/s = 10, 30, 100 
and 300 GeV for GBW (left) and KST (right) parameterizations. 



(ex rf.) for the dipole cross section is used. For the coefficient in front of r\ we employ the 
expression obtained by the first term of Taylor expansion of Eq. (0): 



«»,2 " ■ 



a [s 



°qq(r T ,s) = ~YTX- r T- ( 12 ) 



rg(s) 

2.2 Charmonium wave functions 

The spatial part of the cc pair wave function satisfying the Schrodinger equation 

A 



+ V(r) V nlm (r) = ^mm(r) (13) 



is represented in the form 

W) = *rd(r)-Y lm (9,<p) , (14) 

where r is 3-dimensional cc separation, ty n i(r) and Yi m (9, ip) are the radial and orbital parts 
of the wave function. The equation for radial \l/(r) is solved with the help of the program 
13fl . The following four potentials V(r) have been used (see Fig. |3|): 
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"COR" : Cornell potential p 



V(r) = -* + -= (15) 



r er 



with k = 0.52, a = 2.34 GeV" 1 and m c = 1.84 GeV. 

"BT": Potential suggested by Buchmiiller and Tye |14| with m c = 1.48 GeV. It has a 



similar structure as the Cornell potential: linear string potential at large separations 
and Coulomb shape at short distances with some refinements, however. 



• "LOG" : Logarithmic potential (IB 

V{r) = -0.6635 GeV + (0.733 GeV) log(r • 1 GeV) (16) 
with m c — 1.5 GeV. 

• "POW": Power-law potential [0 

V(r) = -8.064 GeV + (6.898 GeV)(r • 1 GeV) ai (17) 

with m c = 1.8 GeV. 

The shapes of the four potentials is displayed in Fig. fj and differ from each other only at 
large r (> lfm) and very small r (< 0.05 fm) separations. Note, however, that COR and 
POW use m c « 1.8 GeV, while BT and LOG use m c « 1.5 GeV for the mass of the charmed 
quark. This difference will have significant consequences. 

The results of calculations for the radial part ^ n i(r) of the IS and 2S states are depicted in 
Fig. |j. For the ground state all the potentials provide a very similar behavior for r > 0.3 fm, 
while for small r the predictions are differ by up to 30%. The peculiar property of the 2S 
state wave function is the node at r ~ 0.4 fm which causes strong cancelations in the matrix 
elements Eq. ([[]) and as a result, a suppression of photoproduction of ip' relative to Jfip 
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Figure 3: Shapes of the potentials V(r) for the four parameterizations employed 
in this paper. The curves for COR, LOG and POW are normalized at r — 1 fm 
to the value of BT potential. 

2.3 Light-cone wave functions for the bound states 

As has been mentioned, the lowest Fock component |cc) in the infinite momentum frame 
is not related by simple Lorentz boost to the wave function of charmonium in the rest 
frame. This makes the problem of building the light-cone wave function for the lowest 
\cc) component difficult, no unambiguous solution is yet known. There are only recipes 
in the literature, a simple one widely used [JHJ, is the following. One applies a Fourier 
transformation from coordinate to momentum space to the known spatial part of the non- 
relativistic wave function (|14]) , \I/(r) => *&(p), which can be written as a function of the 
effective mass of the cc, M 2 = 4(p 2 + m 2 .), expressed in terms of light-cone variables 



M 2 (a,p T ) 



p 2 T + m 2 c 



(18) 
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Figure 4: The radial part of the wave function \l/„;(r) for the IS and 2S states 
calculated with four different potentials (see text). 



In order to change integration variable pl to the light-cone variable a one relates them via 
M, namely pi = (a — l/2)M(pr, a). In this way the cc wave function acquires a kinematical 
factor 

(V , m 2N3/4 

*(P) ^ ^2 T + Jy /2 ■ = ^KPT) • (19) 



This procedure is used in |20| and the result is applied to calculation of the amplitudes 
(HI). The result is discouraging, since the tp' to J/ip ratio of the photoproduction cross sections 
are far too low in comparison with data. However, the oversimplified dipole cross section 
cr gg( r T) oc r\ has been used, and what is even more essential, the important ingredient of 
Lorentz transformations, the Melosh spin rotation, has been left out. The spin transforma- 
tion has also been left out in the recent publication |JT] which repeats the calculations of 



f20|| with a more realistic dipole cross section which levels off at large separations. This leads 
to suppression of the node-effect (less cancelation) and enhancement of photoproduction. 
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Nevertheless, the calculated ip' to J/0 ratio is smaller than the data by a factor of two. 

The 2-dimensional spinors Xc and Xc describing c and c respectively in the infinite mo- 
mentum frame are known to be related via the Melosh rotation [^2|, |19| to the spinors Xc 
and Xc i n the rest frame: 

X c = R(a,p T )Xc 

X- c = R(l-a,-p T ) X -c, (20) 

where the matrix R(a,px) has the form: 

s m c + aM-i[ax n]p T 

R{a 1 p T ) = , • (21) 

^J(m c + aM) 2 + p T 

Since the potentials we use in section [2.2| contain no spin-orbit term, the cc pair is in 
S'-wave. In this case spatial and spin dependences in the wave function factorize and we 
arrive at the following light cone wave function of the cc in the infinite momentum frame 

^\a,p T ) = U^(a,p T ) ■ ^(a,p T ) , (22) 

where 

U^\a,p T ) = xf R f (a,p T ) B ■ e> a y R*(l - a, -p T ) a~ l %% (23) 

and Xc is defined in (||). 

Note that the wave function (|22|) is different from one used in |2^, [24], |25] where it was 
assumed that the vertex ip — > cc has the structure ip^u^^u like the for the photon 7* — > cc. 
The rest frame wave function corresponding to such a vertex contains 5* wave and D wave. 
The weight of the latter is dictated by the structure of the vertex and cannot be justified by 
any reasonable nonrelativistic potential model for the cc interaction. 

Now we can determine the light-cone wave function in the mixed longitudinal momentum 
- transverse coordinate representation: 

^\a,r T ) = -Ly d 2 p T c-^ ^\a,p T ) . (24) 
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The spatial component <J>^(a,r T ) of Eq. fll9|) in mixed representation (|24D is plotted as 
a function of tt and a in Fig. |5] for J/ip(lS) and ip'(2S) states. While the IS* wave function 
depends monotonically on Tt and smoothly vanishes at small a, the wave function of the 2S 
state demonstrates a nontrivial behavior: the node disappears for small a. 




Figure 5: Three-dimensional plot for the light-cone wave functions for J/i()(lS) 
and i/j'(2S) in the mixed a — f t representation for the BT potential |T^/ . 



3 Calculations and comparison with data 
3.1 The final expressions 

Having the light-cone wave function of charmonium in momentum representation Eq. ( p2[ ) 
it is more convenient to switch to an integration over in the matrix element Eq. ([!]): 

l 

M t ,l( s iQ 2 ) = j da J d 2 p T %(aJ T )^ T A a JT,s,Q 2 ) , (25) 
o 
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where 



1 



Z T . L (a,p T ,s,Q 2 ) = ^-J2 uM (^Pt) / d 2 r T e ip ^ a(r T ,s)^t\a,f T ) 



(26) 

If the dipole cross section depends on tt like cr (l — e~' r T/ r oj (see Eqs. ([5]) and (0)), then 
S(a, pt, s, Q 2 ) which includes the effects of spin rotation, can be expressed as follows: 



S T (a,p T ,s,Q 2 ) = — 



rriT 

2 



2p\ a(l — a) 



2 p' T 



m q rriT r$ 



1 + 



itlt + rriL 
rriT (1 — 2a) 



j: T (a,p T ,s,Q 2 ) 
df: t (a,p T ,s,Q 2 ) 



m T + rriL 



dp\ 



T, L (a,p T ,s,Q 2 



m c (rriT + mi) 
where = m 2 + p\, m 2 L = Am 2 a(l — a) and 

(T (s x 



m z + m T m L ~ 2 
T, L (a,pT,s,Q ) 



(27) 
(28) 



Y, TiL (a,p T ,s,Q 2 ) 
t, t (a,p T ,s,Q 2 ) 

<5> T (a,r T ,Q 2 ) 
$ L (a,r T ,Q 2 ) 



2tt 
2tt 



dVe 1 ^ $ T ,L(«,r- T ,Q 2 )- 
rf 2 r T e^ T $ T (a,r T ,Q 2 ) -e 



1 - e~ r T 



-r|/rg(s) 



rf/rg(s) 



1 2a, 



7T V 3 

2 / C^p.m. 



m q K (er T ) 



7i v 3 

The photoproduction cross section is given by 



Qa(l - a) K (er T ) . 



(29) 
(30) 

(31) 
(32) 



a. 



7*p— j^P 



|.M T (s,Q 2 )| 2 + g |.M L ( g ) Q 

167T.B 



2\|2 



(33) 



where e is the photon polarization (for HI data (e) = 0.99); B is the slope parameter in 
reaction j*p — > ipp. We use the experimental value J29| B = 4.73 GeV~ 2 . JAt,l includes 
also the correction for the real part of the amplitude: 

< 7T d\nM T As,Q 2 )\ 



Mt,l(s,Q) = M t ,l(s,Q) I"* 



(34) 



2 5 1ns J 

where we apply the well known derivative analyticity relation between the real and imaginary 
parts of the forward elastic amplitude ||26|| . The correction from the real part is not small 
since the cross section of charmonium photoproduction is a rather steep function of energy 
(see below). 
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3.2 s and Q 2 dependence of cr{^*p — > J/ipp) 

Now we are in a position now to calculate the cross section of charmonium photoproduction 
using Eq. (^). The results for J/ip are compared with the data in Fig. ^. Calculations are 
performed with GBW and KST parameterizations for the dipole cross section and for wave 
functions of the J/ip calculated from BT, LOG, COR and POW potentials. One observes 

• There are no major differences for the results using the GBW and KST parameteriza- 
tions. 

• The use of different potentials to generate the wave functions of the J/ip leads to two 
distinctly different behaviors. The potentials labeled BT and LOG (see sect. |2.2j) 



describe the data very well, while the potentials COR and LOG underestimate them 
by a factor of two. The different behavior has been traced to the following origin: BT 
and LOG use m c « 1.5 GeV, but COR and POW m c « 1.8 GeV. While the bound 
state wave functions of J/ip are little affected by this difference (see Fig. the photon 
wave function Eq. (|3]) depends sensitively on m c via the argument Eq. fl5|) of the K 
function. 

We compare our calculations also with data for the Q 2 dependence of the cross section. 
The data are plotted in Fig. [7| at cm. energy y/s = 90 GeV as a function of Q 2 + M 2 ^, since 
in this form both, data and calculations display an approximate power law dependence. Such 
a dependence on Q 2 + M 2 ^ is suggested by the variable e 2 in Eq. (||), which for a = 1/2 takes 
the value Q 2 + (2m c ) 2 . It may be considered as an indication that a = 1/2 is a reasonable 
approximation for the nonrelativistic charmonium wave function. 

Our results are depicted for BT and COR potentials and using GBW and KST cross 
sections. Agreement with the calculations based on BT potential is again quite good, while 
the COR potential grossly underestimate the data at small Q 2 . Although the GBW and KST 
dipole cross sections lead to nearly the same cross sections for real photoproduction, their 
predictions at high Q 2 are different by a factor 2 — 3. Supposedly the GBW parameterization 
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Figure 6: Integrated cross section for elastic photoproduction jp — > J/ipp with 
real photons (Q 2 = 0) as a function of the energy calculated with GBW and 
KST dipole cross sections and for four potentials to generate Jfifj wave functions. 
Experimental data points from the HI E401 gfljj, E516 g3jj and ZEUS fft]/ 
experiments. 



should be more trustable at Q 2 ^> M|. 

3.3 Importance of spin effects for the ip' to J/i/j ratio 

It turns out that the effects of spin rotation have a gross impact on the cross section of 
elastic photoproduction •yp — > J/ip(ip)p . To demonstrate these effects we present the results 
of our calculations at ^fs = 90 GeV in Table |l|. The upper half of the table shows the pho- 
toproduction cross sections for J/ip for different parameterizations of the dipole cross section 
(GBW, KST, "rf,") and potentials (BT, COR, LOG, POW). The numbers in parenthesis 
show what the cross section would be, if the spin rotation effects were neglected. We see 
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Figure 7: Integrated cross section for elastic photo production as a function 
of the photon virtuality Q 2 + Mj/^ at energy ^Js = 90 GeV. Solid and dashed 
curves are calculated with GBW and KST dipole cross sections, while thick and 
thin curves correspond to BT and COR potentials, respectively. Results obtained 
with LOG and POW potentials are very close to that curves (LOG similar to BT 
and POW to COR, see also Fig. §). 



that these effects add 30-40% to the J/ifj photoproduction cross section. 

The spin rotation effects turn out to have a much more dramatic impact on ijj' increasing 
the photoproduction cross section by a factor 2-3. This is visible in the lower half of the 
table which shows the ratio R = a{jp') / a(J/ip) of photoproduction cross sections, where the 
number in parenthesis correspond to no spin rotation effects included. This spin effects 
explain the large values of the ratio R observed experimentally. Our results for R are about 
twice as large as evaluated in |^TJ and even more than in [E21J . 

The ratio of ip' to J/0 photoproduction cross sections is depicted as function of cm. 
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Table 1: The photoproduction 7p — > J/tpp cross-section a(J/ip) in nb and the 
ratio R = a{ip') / a(J/ip) for the four different types of potentials (BT, LOG, 
COR, ROW) and the three parameterizations (GBW, KST, r\) for the dipole 
cross section a(rx, s) at s/s = 90 GeV. The values in parentheses correspond to 
the case when the spin rotation is neglected. 





BT 


LOG 


COR 


POW 


[0] GBW 
KST 

r\ 


52.01 (37.77) 
49.96 (35.87) 
66.67 (47.00) 


50.78 (36.63) 
48.49 (34.57) 
64.07 (44.86) 


23.13 (17.07) 
21.05 (15.42) 
25.81 (18.71) 


24.94 (18.64) 
22.83 (16.92) 
28.23 (20.66) 


R GBW 
KST 

r\ 


0.147 (0.075) 
0.147 (0.068) 
0.101 (0.034) 


0.117 (0.060) 
0.118 (0.054) 
0.081 (0.027) 


0.168 (0.099) 
0.178 (0.099) 
0.144 (0.070) 


0.144 (0.085) 
0.152 (0.084) 
0.121 (0.058) 



energy in Fig. ^ and as a function of Q 2 in Fig. |9| for all four potentials and for the parame- 
terizations of the dipole cross sections GBW and KST. 

Our calculations agree with available data, but error bars are too large to provide a more 
precise test for the theory. Remarkably, the ratio R{s) rises with energy. This result is in 
variance with the naive expectation based on the larger size of the ip' and on the usual rule: 
the smaller the size of the qq dipole, the steeper its energy dependence. There is, however, 
no contradiction, since this is another manifestation of the node in the wave function of if)'. 
Indeed, as function of energy mostly the short distance part of the dipole cross section a q q{rT) 
rises. It enhances the positive contribution for distances shorter than the node position in 
the wave function. Therefore, with increasing energy the cancelation in the amplitude of 
if)' production is reduced. This effect leads to a steeper energy dependence of ip' production 
compared to J/ip. The effect is stronger for GBW than KST parameterizations, since the 
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Figure 8: The ratio of ip' to J/ip photoproduction cross sections as a function 
of cm. energy calculated for all four potentials with with GBW and KST pa- 
rameterizations for the dipole cross section. Experimental data points from the 

EMC " 



SLAC JM, NAM Ml, E401 



NMC J33/ and HI |2f experiments. 



GBW cross section does not rise with energy at all at large separations. Note that this 
situation is specific for photoproduction because the nodeless wave function of the photon is 
projected to the sign changing wave function of ip'. This should not happen in the case of 
elastic J/ip(^')-p scattering (see below). 

Similarly of the node effect leads to a rising Q 2 -dependence of the ip' to J/ip ratio in the 
photoproduction cross sections. Our calculations are compared with available data in Fig. |9] 
for the GBW and KST parameterizations respectively. 
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Figure 9: The ratio of ip 1 to J/ip photoproduction cross sections as a function 
of the photon virtuality Q 2 at energy y/s = 90 GeV all four potentials with with 
GBW and KST parameterizations for the dipole cross section. Experimental 



data points from the HI experiment /26|/. 



4 Charmonium-nucleon total cross sections 



After the light-cone formalism has been checked with the data for virtual photoproduction 
we are in position to provide reliable predictions for charmonium-nucleon total cross sections. 
The corresponding expressions are given by Eq. (|2|)) (compare with ||). For the GBW and 
KST dipole cross sections, which have the form a (l — e~ r T/ r o^j (see Eqs. (H) and @), a 
summation over spin indexes in (0) gives for the .S-states, 



1 oo 



Jdot Jdp T Jdq T U(a,p T ) U(a, q T ) V(a,p T , q T ) 



ooo 



where 



U(a,p T ) = p T ^a,p T )e-^T/ i [(M 2 l (p T )+ P 2 T )(M 2 (p T )+p 2 T 



-1/2 



(35) 



(36) 
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V(a,pT,q T ) = M 1 (p T )M 1 (q T )M 2 (p T )M 2 (q T ) I (v) 



(37) 



+ 



Afi(pa 



M 1 (p T )M 1 (q T ) + Mi(pr)M 2 (sr)l p T q T h(v) + P%qlh(y) 
m c + rriT\ 



M 2 (p T ) = m c + m T ^ 
1 2 

v = o r 0PTQT 



a 


1 




a 


1 




a 


a 



(38) 

(39) 
(40) 



Here = m|! + $^(«,pt) is defined in (|I9~D; Jo,i,2(^) are Bessel functions of imaginary 
argument. 



The calculated J ftp- and ^//-nucleon total cross sections are plotted in Fig. |10| for for the 
GBW and KST forms of the dipole cross sections and all four types of the charmonium 
potentials. 
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Figure 10: Total J/tpp (thick curves) and ip' p (thin curves) cross sections with 
the GBW and KST parameterizations for the dipole cross section. 
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The corresponding results for %-states are depicted in Fig. O. Here m is the projection 
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Figure 11: Total XP ( m — — thick curves, m — 1 — thin curves) cross sections 
with the KST and GBW parameterizations for the dipole cross section. 



of the orbital momentum which can be zero or one, since this is a P-wave state. From these 
cross sections with definite m, which we denote u^, one can construct the total cross sections 
for the Xc states with different spins and helicities A, 



Xc0 (A = 0): a = i(2of + 

Xd(A = 0) : o = o\- 

Xci(A = ±l): a = + < 

Xc2 (A = 0): a = 1 (<7? + 2<r8 ; , 

Xc2 (A = ±l): a = l - (a 1 x + a x ) ; 

Xc2 (A = ±2): o = o\. (41) 
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Using these relations one can easily derive the cross sections averaged over helicities which 
are equal for all three states Xc 0,1,2- 

The strong dependence of the cross sections for P-wave charmonium states on the pro- 



jection m = 0, 1 of the orbital momentum has been found previously in pOfl . However the 
predicted cross sections at s/s = lOGeV for Xc( m = 0), Xcijn = 1) and ip' are about twice 
as large as ours. We believe that the disagreement originates from the too rough nonper- 
turbative dipole cross section^ used in j|(J which was not well adjusted to data. Even the 



pion-nucleon cross section calculated with Eq. (1) in ||0| overestimates the experimental 
value by factor 1.5. 

Although all four potentials are presented, comparison with photoproduction data in 
Figs. ^| and [7] show that two of them, BT and LOG potentials, are more trustable at least for 
J/0. These two potentials again give very close predictions for J/ip-p total cross sections but 
the deviation from the predictions with the two other potentials, COR and POW, is much 
smaller than in the case of photoproduction. 

Note that the cross sections calculated with the GBW parameterization demonstrate 
a tendency to level off at very high energy, especially for ip', as compared to the KST 
predictions. The reason is obvious: the GBW cross sections approach the universal limit 
&max — o"o = 23.03 mb. This cannot be true, and the KST parameterization is more reliable 
than GBW at high energies where the gluon cloud surrounding the cc pair becomes nearly 
as big as light hadrons. 

According to Figs. [TI] and [ll| for the KST parameterization the total cross sections of 
charmonia are nearly straight lines as function of \fs in a double logarithmic representation, 
though with significantly different slopes for the different states. Therefore a parameteriza- 
tion in the form 

= *t'(4-Y » (42) 



1- we are thankful to Lars Gerland who provided us with the expression for the dipole cross section used 
in Bell. 
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seems appropriate, at least within a restricted energy interval. We use the data shown in 
Figs. and |ll] for the KST parameterization of a q q and for the BT and LOG potentials and 



fit the them by the form (|42j) with sq = 1000 GeV. The two values from the BT and LOG 
potentials have been averaged and their half difference gives the error estimation. Table ||] 
shows values for Oq and A averaged over the energy interval 10 GeV < \fs < 300 GeV, 
and the bound state sizes (rfi). As expected erg rises monotonically with the size of the 

Table 2: Averaged sizes (r^) for charmonia bound states together with cr and 
A in the parameterization ( p^ j for the J/ip-, ip'- and x-P r °t° n cross sections. 
Estimation of the errors is given in the text. 





(r 2 ) [fm 2 ] 


at [mb] 


A 




J/i> 




0.117 ±0.003 


5.59 ±0.13 


0.212 ±0.001 


X 


[m = 


0) 


0.181 ±0.004 


7.17 ±0.07 


0.195 ±0.001 


X 


[m = 


1) 


0.362 ±0.007 


13.17 ±0.16 


0.164 ±0.002 








0.517 ±0.034 


16.63 ±0.59 


0.139 ±0.005 



charmonium state, and the cross section for if)' N is about three times larger that for J/ip. 
This deviates from the r 2 scaling, since the mean value (r 2 ) is 4 times larger for ip' than for 
J/ijj. The exponent A which governs the energy dependence decreases monotonically with 
the size of the charmonium state, demonstrating the usual correlation between the dipole size 
and the steepness of energy dependence. The values of A are larger than in soft interactions 
of light hadrons (~ 0.08), but smaller than values reached in DIS at high Q 2 . 

Our results at \fs = 10 GeV (the mean energy of charmonia produced in the NA38 /NA50 
experiments at SPS, CERN), 

aS{Vs = 10 GeV) = 3.56 ± 0.08 mb , (43) 
vfotiVs = 10 GeV) = 12.19 ± 0.61mb , (44) 
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well agree with the cross sections extracted in from photoproduction data employing 
the two-channel approximation, 2.8 ± 0.12 mb < (Ttlt'is/s = lOGeV) < 4.1 ± 0.15 mb and 
(jfot/o'tot' ~ 3.75 (having poorly controlled accuracy), which shows that the two channel 
approach is a reasonable tool to analyze photoproduction data. 

The cross section Eq. (|4"2"D) with the parameters in Table ^ agrees well with a^f (a/s = 
20 GeV) = 4.4 ± 0.6 mb obtained in the model of the stochastic vacuum jOJ]. 

It worth noting that the results for charmonium-nucleon total cross sections are amazingly 

similar to what one could get without any spin rotation, 

l 

<4<*(?) » Jda J d 2 r T \^(a, r T )f a q q(r T ,s) , (45) 
o 

where $,/,(«, ty) is related by Fourier transformation to Eq. (|19"D , or even performing a 
simplest integration using the nonrelativistic wave functions ( |13|) in the rest frame of the 
charmonium: 

4*{s) « j#r |^(r-)| 2 a g ,(r T , s) . (46) 
The comparison presented in Fig. O for the BT potential shows that fl4"5|) - (^6]) are only 



about 10% below the exact calculation for J/ip, while there is practically no difference 
between the exact and approximate calculations for ip' . 

5 Nuclear suppression of charmonium production 

Production of charmonia off nuclei seem to be a natural source of information about charmonium- 
nucleon cross section since nuclear absorption leads to suppression of the production rate 
measured experimentally. However, one should be cautious applying our results to a cal- 
culation of nuclear attenuation of charmonium. In exclusive photoproduction of charmonia 
the uncertainty principle does not allow to resolve between J/ ip and ip' unless the formation 
time, t f = 2E^/(M%, - Mj^) 0, ||, is shorter than the mean inter-nucleon separation in 



nuclei. Only one experiment j|3| at ~ 20 GeV satisfies this condition. Analyzed with an 
optical model it leads to <jfl^ N = 3.5 ± 0.8 mb in a good agreement with our calculations. 
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Figure 12: Comparison of the results for <jf ^ (s) obtained with the exact ex- 
pression (TJ) (solid curves) and with the approximations ($3() (dashed) and 
(dotted). 



The nuclear photoproduction data ||4j taken at 120 GeV cannot be treated in the same way 
since the formation time, If ~ 10 fm exceeds the nuclear size. In addition, the coherence 
length, l c = 2E^,/(Mj,^ + Q 2 ) || ^IJ, is also long, about 5 fm substantially increasing the 
attenuation path for the produced cc pair. 

In the case of hadroproduction of charmonia off nuclei, the interplay of the formation 
and coherence time effects are as important as in photoproduction. On top of that, the 
situation is complicated by decays of x s an d i>' which substantially feed the yield of J/ip. 
These heavier states, even if their absorption cross sections are known from our calculations, 
are also subject to the effects of formation and coherence lengths. 



In the analysis \ W\ of data from the experiment E866 of pA — > J/ipX collisions at 800 GeV, 
proper attention has been given to coherence and formation time effects with the result 
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(extrapolated to yfs = 10 GeV) 



0", 



"Wp 

eff 



5.0±0.4mb , 



(47) 




10.5 ± 3.6 mb 



(48) 



The effective J/-?/;- nucleoli cross section which is fed by decays of heavier states can be 
estimated as follows, 



where (T) « 0.75 paRa is the mean thickness of a nucleus with radius Ra and the mean 
density pa ~ 0.16 fm~ 3 . 

Eq. (|4~9"D is relevant for Jfijj suppression in nuclear collisions (proton-nucleus and nucleus- 
nucleus). In this reactions the observed J/ip arises from directly produced J/0's with prob- 
ability w\ < 1 and from the other states x> ^ y i a decay after the charmonia have left the 
interaction zone, where W{ is the probability that the state "i" contributes to the finally 
observed J/0. Values for Wi and afjf are given in Table ^ where m — 0, 1 is the projection 
of the orbital momentum of the cc pair on the direction of gluon-gluon collision in Xifl pro- 
duction (xo has a tiny branching to J/ip). It turns out that Xi and X2 with m = cannot be 
produced or are strongly suppressed in gluon fusion due to the selection rules which forbid 
projections ±1 for the total angular momentum (e.g. see |47|), this is why we put W2 = 0. 

We calculate a e ff for tungsten used in the analysis |46| and find for ^fs = 10 GeV 



where the main uncertainty arises from the Wi. This number is in a good accord with 
Eq. (fI7D, while the calculated value for af f Eq. ( |4"4] ) agrees well with (^Sj). 

The coherence effects are quite important even at the energy of the NA38/NA50 exper- 
iments (Ei/, ~ 50 GeV) at CERN, this is why the effective absorption cross section for ip' 
production suggested by the data is about a half of the value we predict. At the energies 
of RHIC and LHC both the coherence and formation times substantially exceed the sizes of 
heavy nuclei, and shadowing becomes the dominant phenomenon. 




(49) 



0. 



."Wp 

eff 



5.8±0.2mb , 



(50) 
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Table 3: Values for the J/ip-, ip'-, x~ an d effective "J/ip" -proton cross sections at 
energy y/s = 10 GeV. Errors are given by averaging on BT and LOG potentials 
for the wave functions. 





Wi 


o~[mb] 








0.52 - 0.6 


3.56 ±0.08 


X 


[m = 


0) 





4.66 ±0.06 


X 


[m = 


1) 


0.32 - 0.4 


9.05 ±0.16 




i>' 




0.08 


12.19 ±0.61 



6 Conclusion and discussions 

In this paper we have proposed a simultaneous treatment of elastic photoproduction o" 7 * p ^^ p (s, Q 2 ) 
of charmonia and total cross sections ufot(s). The ingredients are (i) the factorized light-cone 
expressions ([I]) - (0) for the cross sections; (ii) the perturbative light-cone wave functions for 
the cc component of the 7*; (iii) light cone wave functions for the charmonia bound states, 
and (iv) a phenomenological dipole cross section a qq -(rT, s) for a cc interacting with a proton. 

The dipole cross section rises with energy; the smaller is the transverse qq separation, the 
steeper is the growth. The source of the energy dependence is the expanding cloud of gluons 
surrounding the qq pair. The gluon bremsstrahlung is more intensive for small dipoles. The 
gluon cloud can be treated as a joint contribution of higher Fock states, \qqnG), however, it 
can be also included into the energy dependence of <jq q (rT, s), as we do, and this is the full 
description. Addition of any higher Fock state would be the double counting. 

As function of energy the initial size of the qq source is gradually "forgotten" after multi- 
step radiation, the small cross sections grow steeper and eventually approach the larger ones 
at very high energies. All the cross sections are expected to reach a universal asymptotic 
behavior which saturates the Froissart bound. 
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The effective dipole cross section cr q q(rT, s) is parameterized in a form which satisfies the 
expectations a qq - oc r|> for tt — > (color transparency), but levels off for tt — > oo. Two 
parameterizations for <J qq (rT, s), whose form and parameters have been fitted to describe 
(Ttot( s ) an< i t ne structure function F 2 (x, Q 2 ) are used in our calculations. 

While the description of the photon wave function is quite certain, the light-cone wave 
function of charmonia is rather ambiguous. We have followed the usual recipe in going 
from a nonrelativistic wave function calculated from a Schrodinger equation to a light cone 
form. We have included the Melosh spin rotation which is often neglected and found that 
it is instrumental to obtain agreement, since no parameter is adjustable. In particular, it 
increases the ip' photoproduction cross section by a factor 2-3 and rises the ip' to J / ip ratio 
to the experimental value. 

At the same time, the charmonium-nucleon total cross sections (J/ip, ip', x{ m = 0) an d 
%(m = 1)) turn out to be rather insensitive to the way how the light-cone wave function 
is formed, even applying no Lorentz transformation one arrives at nearly the same results. 
This is why we believe that the predicted charmonium-nucleon cross section are very stable 
against the ambiguities in the light-cone wave function of charmonia. A significant energy 
dependence is predicted which varies from state to state in accord with our expectations. 

We show our predictions for charmonium-nucleon cross sections in a restricted energy 
range 10 GeV < y/s < 300 GeV, but this interval can be largely extended in both directions. 
Since the OZI rule suppresses the leading Reggeons, one can stay with gluonic exchanges 
rather far down to low energies, unless the charmed Reggeon exchanges become important 
||. 
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